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This  work  presents  axial  ion  energy  distribution  measurements  within  the  acceleration 
channel  of  the  600  W  Busek  Co.  Inc.  BHT-HD-600  laboratory  Hall  thruster  derived  from 
laser-induced  fluorescence  measurements  of  the  5d[4]7/2  —  6p[3]5/2  xenon  ion  excited  state 
transition.  Acceleration  channel  centerline  ion  energy  distributions  are  measured  for  three 
closely  related  operating  conditions  which  only  differ  in  the  magnitude  of  the  radial  magnetic 
field  strength.  These  three  operating  conditions  span  a  broad  range  of  discharge  current 
oscillations  strengths.  The  0  to  200  kHz  frequency  domain  is  characterized,  and  the  dominant 
40  kHz  to  50  kHz  frequency  appears  most  likely  to  be  axially  traveling  ionization  waves, 
commonly  known  as  the  breathing  mode  oscillations.  These  oscillations  manifest  themselves 
to  the  laser  induced  fluorescence  diagnostic  as  clearly  broadened  ion  energy  distributions 
when  the  oscillation  strength  is  high.  We  determine  the  spatial  extent  of  the  axial  breathing 
mode  oscillation  nonintrusively.  The  coherence  and  magnitude  of  the  discharge  current 
oscillations  are  inversely  proportional  to  acceleration  channel  radial  magnetic  strength. 


Introduction 

THe  goal  of  this  study  is  to  characterize  the  xenon 
ion  energy  distributions  inside  the  acceleration 
channel  of  a  low  power  xenon  Hall  thruster  using 
laser-induced  fluorescence  (LIF)  and  assess  the  effect 
of  discharge  plasma  oscillations  on  the  distribution 
of  ion  energies.  Although  the  discharge  of  a  Hall 
thruster  appears  to  be  a  quiescent  plasma,  it  can  be 
a  highly  oscillatory  plasma  with  a  variety  of  axial,  az¬ 
imuthal,  and  possibly  radial  modes.  Understanding 
these  oscillations  is  necessary  to  the  application  and 
interpretation  of  the  any  internal  measurements  for 
thruster  plasma  model  validation. 

Ion  velocity  and  energy  distribution  measurements 
are  valuable  for  their  use  in  the  validation  of  device 
models  of  this  and  other  Hall  thrusters  as  they  pro¬ 
vide  a  direct  comparison  of  measurement  to  simula¬ 


tion  results.  Previous  attempts  to  compare  bulk  ion 
velocities  from  experiments  with  results  of  numerical 
simulations  have  encountered  difficulties.  For  exam¬ 
ple,  simulation  efforts  usually  report  mean  velocities. 
Experimental  methods  such  LIF  are  often  noise  lim¬ 
ited  and  are  better  suited  to  determining  the  most 
probable  velocity.  In  the  skewed,  non-symmetric  ve¬ 
locity  distributions  common  to  Hall  thrusters,  the 
most  probable  (i.e.  peak  signal)  and  statistical  mean 
velocities  often  differ.  Ideally,  comparison  of  the 
velocity  distribution  function  determined  from  each 
method  would  minimize  ambiguity.  Due  to  the  rela¬ 
tively  narrow  transition  width  and  wide  velocity  dis¬ 
tribution,  these  fluorescence  traces  have  been  previ¬ 
ously  shown  to  be  representative  of  the  ion  velocity 
distribution  function  and  can  also  be  used  to  calcu¬ 
late  ion  energy  distributions.^ 
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Most  previous  internal  measurements  of  ion  veloc¬ 
ities  and  energies  have  been  carried  out  using  a  slot 
cut  into  the  side  of  a  Hall  thrusters.^’ ^  These  slots 
presumably  affect  the  operation  of  the  thruster.  Even 
if  the  global  effect  is  small,  the  local  effect  may  be  sig¬ 
nificant.  Furthermore,  slicing  slots  into  Hall  thruster 
side  walls  is  not  always  possible  with  expensive,  one 
of  kind,  test  articles. 

Most  modern  Hall  thrusters  have  acceleration 
channels  with  a  maximum  depth  of  1-2  cm.  There¬ 
fore,  it  is  possible  to  align  collection  optics  to  the 
probe  beam  such  that  limited  internal  optical  access 
is  possible  without  modification  of  the  Hall  thruster. 
In  this  work,  the  collection  lens  is  placed  60°  off  the 
plume  axis.  This  minimizes  plume  impact  since  typ¬ 
ically  greater  than  95%  of  a  Hall  thruster  plume  ion 
flux  is  contained  in  a  45°  half  angle.  In  this  way,  it  is 
possible  to  completely  non-intrusively  probe  internal 
ion  acceleration  of  any  modern  Hall  thruster^ 

In  this  study  acceleration  channel  centerline  ion 
energy  distributions  are  measured  for  three  closely 
related  operating  conditions  which  only  differ  in 
the  magnitude  of  the  radial  magnetic  field  strength. 
These  three  operating  conditions  span  a  broad  range 
of  discharge  current  oscillations  strengths  from  mild 
broad  band  oscillations  at  the  highest  radial  magnetic 
field  strength  to  highly  coherent  oscillations  which 
dominate  the  discharge  current  at  the  lowest  mag¬ 
netic  field  strength. 

Experimental  Apparatus 

Xenon  Ion  Spectroscopy 

LIE  is  a  convenient  diagnostic  for  the  investigation 
of  ion  and  atomic  velocities  as  it  does  not  perturb 
the  plasma.  The  LIE  signal  is  a  convolution  of  the 
velocity  distribution,  transition  line  shape,  and  laser 
beam  frequency  profile.  Determination  of  the  velocity 
distribution  from  LIE  data  only  requires  the  decon¬ 
volution  of  the  transition  line  shape  and  laser  beam 
profile  from  the  raw  LIE  signal  trace.  Subsequently, 
it  is  straightforward  to  calculated  the  energy  distri¬ 
bution  from  the  the  velocity  distribution. 

For  the  results  reported  here,  the  5d[4]7/2  —  6p[3]5/2 
electronic  transition  of  Xe  H  at  834.7  nm  is  probed. 
The  isotopic  and  nuclear-spin  effects  contributing 
to  the  hyperfine  structure  of  the  5d[4]Y/2  —  6p[3]5/2 
xenon  ion  transition  produce  a  total  of  19  isotopic 
and  spin  split  components.  The  hyperfine  splitting 
constants  which  characterize  the  variations  in  state 
energies  are  only  known  for  a  limited  set  of  energy 
levels.  Unfortunately,  the  834.7  nm  xenon  ion  tran¬ 
sition  only  has  confirmed  data  on  the  nuclear  spin 
splitting  constants  of  the  6p[3]5/2  upper  state. 
Manzella  first  used  the  5d[4]7/2  —  6p[3]5/2  xenon  ion 


transition  at  834.7  nm  to  make  velocity  measure¬ 
ments  in  a  Hall  thruster  plume. ^  A  convenient  fea¬ 
ture  of  this  transition  is  the  presence  of  a  relatively 
strong  line  originating  from  the  same  upper  state 
(6s [2] 3/2  —  6p[3]5/2  transition  at  541.9  nm,^  which 
allows  for  non-resonant  fluorescence  collection).  Ion 
velocity  is  simply  determined  by  measurement  of  the 
Doppler  shift  of  the  absorbing  ions.^^ 

Previous  measurements  and  analysis  have  shown 
that  deconvolution  is  not  strictly  required  to  estimate 
xenon  ion  velocity  and  energy  distributions  from  the 
raw  LIE  data  in  this  plasma  discharge  for  this  par¬ 
ticular  xenon  transition.^  The  transition  is  relatively 
narrow  (approximately  600  MHz)  and  the  velocity 
and  energy  distributions  in  the  vicinity  of  the  exit 
plane  is  sufficiently  broad  that  the  fluorescence  trace 
does  not  require  deconvolution  of  the  transition  line 
shape  to  produce  adequate  estimates  of  the  distribu¬ 
tions.  Not  performing  the  deconvolution  further  in 
the  plume  (e.g.  beyond  the  cathode  plane)  may  in¬ 
troduce  uncertainties  estimated  to  be  less  than  20%. 
Zeeman  splitting  is  neglected  in  this  analysis  due  to 
the  broad  velocity  distributions. 

Test  Facility 

The  LIE  measurements  were  performed  in  Chamber 
6  of  the  Air  Force  Research  Laboratory  (AFRL)  Elec¬ 
tric  Propulsion  Laboratory  at  Edwards  AFB,  CA. 
Chamber  6  is  a  non-magnetic  stainless  steel  cham¬ 
ber  with  a  1.8  m  diameter  and  3  m  length.  It  has 
a  measured  pumping  speed  of  32,000  1/s  on  xenon. 
Pumping  is  provided  by  four  single  stage  cryogenic 
panels  (single  stage  cold  heads  at  25  K)  and  one 
50  cm  two  stage  cryogenic  pump  (12  K).  Chamber 
pressure  during  thruster  operation  is  approximately 
1.5  X  10“^  Torr,  corrected  for  xenon. 

The  thruster  is  mounted  on  a  three  axis  orthog¬ 
onal  computer  controlled  translation  system.  Fig¬ 
ure  I  shows  the  Hall  thruster  and  optics  mounted 
within  the  vacuum  chamber  as  well  as  the  LIE  appa¬ 
ratus.  The  laser  is  a  tunable  diode  laser.  It  is  capa¬ 
ble  of  tuning  approximately  ±50  GHz  about  a  center 
wavelength  of  834.7  nm.  The  6  mW  beam  is  passed 
through  a  Faraday  isolator  to  eliminate  feedback  to 
the  laser.  The  laser  beam  then  passes  through  several 
beam  pick-offs  until  it  is  focused  by  a  lens  and  enters 
the  vacuum  chamber  through  a  window.  The  probe 
beam  is  chopped  at  a  frequency  by  an  optical  chopper 
(Ch2  at  2.8  kHz)  for  phase  sensitive  detection  of  the 
fluorescence  signal. 

The  two  wedge  beam  pick-offs  (BS)  shown  in  Fig.  I 
provide  portions  of  the  beam  for  diagnostic  purposes. 
The  first  beam  pick-off  directs  a  beam  to  a  photo¬ 
diode  detector  (DI)  used  to  provide  constant  power 
feedback  to  the  laser.  The  second  beam  is  divided 
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Figure  1.  Diagram  of  the  laser  optical  train  and  collection 
optics. 


into  two  equal  components  by  a  50-50  cube  beam 
splitter.  The  first  component  is  directed  to  a  wave¬ 
length  meter  used  to  monitor  absolute  wavelength. 
The  second  component  is  sent  through  an  optical 
chopper  (Chi  at  1.4  kHz)  and  through  a  low  pressure 
xenon  hollow  cathode  discharge  lamp.  The  lamp  pro¬ 
vides  a  stationary  absorption  reference  for  the  deter¬ 
mination  of  the  Doppler  shift.  Unfortunately,  there  is 
no  detectable  population  of  the  ionic  xenon  5d[4]7/2 
state.  However,  there  is  a  nearby  (estimated  to  be 
18.1  GHz  distant)  neutral  xenon  65'[l/2]i  —  6p'[3/2]2 
transition  at  834.68  nm.^^’^^  The  second  pick-off 
sends  a  beam  to  a  300  MHz  free  spectral  range  Fabry- 
Perot  etalon  (F-P)  that  provides  high  resolution  fre- 


Table  1.  Nominal  BHT-HD-600  Hall  thruster  operating 
conditions. 


Anode  Flow 

2.453  mg/s 

Cathode  Flow 

197  /ig/s 

Anode  Potential 

300  V 

Anode  Current 

2.16  A 

Keeper  Current 

0.5  A 

Heater  Current 

3.0  A 

quency  monitoring  of  the  wavelength  interval  swept 
during  a  laser  scan. 

The  fluorescence  collection  optics  are  also  shown  in 
Fig.  1.  The  fluorescence  is  collected  by  a  75  mm  di¬ 
ameter,  300  mm  focal  length  lens  within  the  chamber 
and  oriented  60°  from  the  probe  beam  axis.  The  col¬ 
limated  fluorescence  signal  is  directed  through  a  win¬ 
dow  in  the  chamber  side  wall  to  a  similar  lens  that 
focuses  the  collected  fluorescence  onto  the  entrance 
slit  of  the  125  mm  focal  length  monochrometer  with 
a  photomultiplier  tube  (PMT).  Due  to  the  1:1  magni¬ 
fication  of  the  collection  optics,  the  spatial  resolution 
of  the  measurements  is  determined  by  the  geometry 
of  the  entrance  slit  0.7  mm  width  and  1.5  mm  height 
as  well  as  the  sub-mm  diameter  of  the  probe  beam. 
This  apparatus  allows  for  limited  probing  of  the  in¬ 
terior  acceleration  channel  of  Hall  thrusters  with  rel¬ 
atively  shallow  acceleration  channels.  Measurements 
show  that  this  combination  of  apparatus  and  laser 
power  are  well  within  the  linear  fluorescence  regime. 

Hall  Effect  Thruster 

The  Hall  thruster  used  in  this  study  is  the  rectangular 
600  W  Busek  Company  BHT-HD-600  Hall  thruster 
with  a  Busek  3.2  mm  hollow  cathode.  A  photo¬ 
graph  of  the  rectangular  BHT-HD-600  Hall  thruster 
is  shown  in  Fig.  2.  This  thruster  has  a  conventional 
five  magnetic  core  (one  inner,  four  outer)  magnetic 
circuit.  The  acceleration  channel  of  the  thruster  has 
a  32  mm  outer  radius  and  a  24  mm  inner  radius.  The 
acceleration  channel  has  a  depth  of  approximately 
10  mm  between  the  geometrical  exit  plane  and  the 
furthest  forward  extent  of  the  anode.  The  thruster 
has  been  extensively  characterized  to  have  a  thrust 
of  39  mN  with  a  specific  impulse  of  1,530  s  yielding 
an  efficiency  of  approximately  50%  at  the  nominal 
conditions  specified  in  Table  1. 

Figure  3  shows  a  cut-a-way  view  of  the  near  field 
geometry  of  the  Busek  BHT-HD-600  Hall  thruster. 
The  locations  of  the  central  magnetic  poles  and  edges 
of  the  acceleration  channel  are  indicated  as  is  the  po¬ 
sition  of  the  cathode.  The  Cartesian  coordinate  sys¬ 
tem  and  origin  used  in  these  measurements  are  also 
shown  in  Fig.  3,  with  positive  X  axis  going  into  the 
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Figure  2.  Photograph  of  the  Busek  BHT-HD-600  Hall  ef¬ 
fect  thruster. 


Figure  3.  Cross-section  of  the  BHT-HD-600  Hall  thruster 
showing  coordinate  system  of  LIF  measurements. 


page.  The  origin  is  at  the  center  of  the  thruster  face 
due  to  the  ease  and  repeatability  with  which  this  posi¬ 
tion  could  be  located.  All  measurements  presented  in 
this  work  will  lie  on  the  Z  =  0  mm  and  Z  =  28  mm 
line.  This  corresponds  to  the  acceleration  channel 
centerline  opposite  the  cathode  in  the  Y-Z  plane. 

Results  and  Discussion 

Magnetic  Field  Variation  Effects 

The  magnetic  field  of  the  Hall  thruster  may  be  var¬ 
ied  by  changing  the  applied  current  to  either  of  the 
four  series  connected  outer  magnetic  cores,  and,  or 
the  central  core.  In  addition  to  the  nominal  case, 
two  additional  cases  are  investigated  in  this  effort. 
These  cases  represent  geometrically  similar  radial 
magnetic  field  profiles  with  decreases  in  the  magni¬ 
tude  to  85%  and  50%  relative  to  the  nominal  case. 
Figure  4  shows  the  calculated  shape  of  the  radial 
magnetic  field  along  the  acceleration  channel  center 
(Y  =  28  mm,  X  =  0  mm).  Although  the  two  off  nom¬ 
inal  cases  are  15%  and  50%  less  in  magnitude  than 
nominal,  the  profiles  are  geometrically  similar.  The 
two  off-nominal  cases  show  a  maximum  -5%  deviation 
from  the  nominal  radial  B  field,  primarily  between 
Z  =  -4  mm  and  -8  mm. 

In  the  off-nominal  cases,  varying  the  magnetic  field 
does  not  change  the  average  anode  current  apprecia¬ 
bly.  The  change  from  nominal  is  approximately  2%. 


Figure  4.  Relative  radial  magnetic  field  along  the  acceler¬ 
ation  channel  center. 


The  effect  of  magnetic  field  variation  on  performance 
has  yet  to  be  characterized  as  thrust  measurements 
for  the  two  off  nominal  conditions  studied  in  this  work 
are  not  yet  available. 

Other  measures  of  the  plume  such  as  divergence 
measurements  from  Faraday  probe  measurements 
show  differences,  but  these  are  appear  to  be  due  to 
chamber  effects  as  the  measurement  were  taken  at  ei¬ 
ther  Busek,  or  AFRL.  The  largest  published  set  of 
plume  data,  by  Ekholm,  et  al,  looks  extensively  at 
the  effect  of  plume  divergence  due  to  variations  in 
anode  flow  rates  and  discharge  potentials,  but  does 
not  examine  the  effect  of  changes  to  the  magnetic 
field  strength.  It  is  anticipated  that  a  future  effort 
will  measure  the  plume  divergence  differences  of  these 
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Figure  5.  Nominal  anode  current  time  and  frequency  do¬ 
main  behavior. 


-lOOjLvs  -50  0  50  100 

Time 


Figure  6.  Anode  current  time  and  frequency  domain  be¬ 
havior  for  85%  nominal  magnetic  field  case. 


three  cases. 

One  effect  that  is  clearly  identified  and  the  primary  ^ 
reason  for  the  study  of  these  three  particular  cases  g 

(despite  the  limited  performance  and  plume  data)  u 

is  the  distinct  anode  current  fluctuation  behavior  of  "§ 
each  case.  Figure  5  shows  the  anode  current  time  < 
and  frequency  domain  behavior  of  the  nominal  mag¬ 
netic  field  case.  The  anode  current  is  relatively  qui¬ 
escent  with  a  mean  of  2.16  A  and  standard  deviation 
of  77  mA  (3.6%). 

Reducing  the  magnetic  field  by  15%  increases  the 
current  fluctuations  to  approximately  5%  as  shown 
in  Fig.  6.  Here  the  mean  anode  current  is  2.21  A 
with  a  standard  deviation  of  112  mA.  These  two  rel¬ 
atively  quiescent  cases  can  be  contrasted  with  the 
50%  nominal  magentic  field  case  in  Fig.  7  which  has 
an  anode  current  fluctuation  standard  deviation  of 
over  34%  of  the  2.10  A  mean.  Although  the  shape 
of  the  magnetic  field  does  not  change  appreciably, 
the  anode  current  oscillations  are  strongly  influenced. 

The  characteristics  of  these  anode  current  fluctua¬ 
tions  are  quantitatively  illustrated  in  the  frequency 
domain  plots  also  shown  in  Figs.  5,  6,  and  7.  As 
the  magnitude  of  the  radial  magnetic  field  shown  in 
Fig.  4  decreases,  the  magnitude  of  the  oscillations  in¬ 
creases.  Interestingly,  as  the  oscillations  increase  in 
magnitude,  they  also  increase  in  coherence.  A  sin¬ 
gle  frequency  (40  kHz)  dominates  the  50%  magnetic 
field  case,  where  as  the  frequencies  of  oscillation  for 
the  nominal  case  are  more  broadly  dispersed. 

It  has  been  hypothesized  that  the  turbulent  plasma 
oscillations  in  the  plasma  are  responsible  for  either 
all  or  at  least  a  portion  of  cross  field  electron  diffu- 
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Figure  7.  Anode  current  time  and  frequency  domain  be¬ 
havior  50%  nominal  magnetic  field  case.. 


sion,  also  commonly  referred  to  as  anomalous  electron 
transport There  are  a  number  of  turbulent  modes 
that  could  account  for  the  time  and  frequency  domain 
behavior  seen  in  Figs.  5,  6,  and  7d^A6  xhese  include 
both  azimuthal  and  axial  ionization  waves  that  may 
in  certain  cases  strongly  influence  the  locations  of  ion 
creation  and  subsequent  acceleration  within  the  Hall 
thruster.  This  waves  are  known  to  occur  at  frequen¬ 
cies  ranging  between  10  kHZ  and  60  kHz,  similar  to 
those  shown  above. 
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Figure  8.  Energy  distribution  functions  between  the  anode  and  geometrical  thruster  exit  plane  along  the  acceleration 
channel  centerline. 
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Energy  distribution  functions  outside  the  geometrical  thruster  exit  plane  along  the  acceleration  channel 


It  is  assumed  in  this  work  that  the  oscillations  seen 
in  this  work  are  primarily  axial  and  correspond  to  the 
so-called  breathing  mode.  This  assumption  appears  to 
be  consistent  with  previous  results  which  have  shown 
the  azimuthal  wave  to  dominate  at  relatively  low  dis¬ 
charge  potentials  (less  than  200  V)  and  dissipate  at 
higher  values. Since  the  anode  oscillations  are  in 
the  correct  frequency  and  operating  regime,  it  ap¬ 
pears  reasonable  to  consider  the  oscillations  to  be 
indicative  of  axially  traveling  ionization  and  poten¬ 
tial  waves  which  strongly  influence  the  local  plasma 
density  and  conductivity. 

LIF  Measurements 

If  the  plasma  potential  and  ionization  locations  are 
axially  traversing  portions  of  the  acceleration  channel 


with  frequencies  on  the  order  of  10 ’s  of  kHz,  the  rela¬ 
tively  long  time  averaged  velocity  distributions  (typ¬ 
ically  4-7  minutes)  will  capture  approximately  of  10^ 
of  these  events.  As  such,  the  resulting  LIF  trace  can 
be  used  to  nonintrusively  characterize  the  severity  of 
the  oscillations  as  well  as  their  physical  extent.  The 
three  cases  for  which  the  oscillatory  behavior  is  char¬ 
acterized  in  Figs.  5,  6,  and  7  represent  an  ideal  test 
case  for  LIF  investigation  of  internal  Hall  thruster 
oscillations. 

Figures  8  and  9  present  the  energy  distribu¬ 
tion  functions  ranging  from  just  above  the  anode 

(Z  =  -8  mm),  through  the  geometrical  exit  plane 

(Z  =  0  mm),  to  10  mm  beyond  the  exit  plane 

(Z  =  10  mm).  All  these  measurements  were  taken 

at  the  acceleration  channel  center  (Z  =  28  mm)  on 
the  X-Y  plane  (X  =  0  mm).  Since  the  plasma  oscil- 
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lations  are  the  result  of  plasma  potential  variations, 
we  have  chosen  to  present  the  LIF  data  in  terms  of 
the  ion  energy  distribution  function  (EDF).  Scaling 
for  all  the  traces  is  by  area  normalization  to  unity. 

In  Figs.  8  and  9,  the  blue  traces  correspond  to  the 
nominal  low  oscillation  case,  the  green  traces  corre¬ 
spond  to  the  80%  radial  magnetic  field  strength  case, 
and  the  red  traces  correspond  to  the  highly  oscilla¬ 
tory  50%  radial  magnetic  field  strength  case.  Exam¬ 
ining  Eigs.  8  and  9,  it  is  immediately  obvious  that 
the  oscillations  have  a  strong  influence  on  the  energy 
distribution.  The  nominal  case  and  the  85%  nominal 
magnetic  field  strength  case  are  are  most  similar.  A 
portion  of  the  differences  in  relative  peak  hight  can 
be  attributed  to  uncertainties  arising  from  the  selec¬ 
tion  of  baselines  during  the  normalization  procedure 
due  to  the  noise  in  the  signal.  Uniformly,  the  low 
magnetic  field  case  (50%  nominal)  energy  distribu¬ 
tions  are  two  to  three  times  broader  than  the  nomi¬ 
nal  and  85%  cases.  This  is  particularly  pronounced 
in  the  interior  of  the  acceleration  channel.  Eor  ex¬ 
ample  at  Z  =  -3  mm,  the  energy  distribution  of  the 
low  magnetic  field  strength  case  varies  between  0  eV 
and  180  eV.  This  is  approximately  3  times  the  width 
of  the  profiles  of  the  other  two  higher  magnetic  field 
cases. 

The  energy  profiles  for  each  case  are  shown  in 
Eig.  10.  Here  for  each  case  we  have  the  median  veloc¬ 
ity  and  width  of  the  profile  calculated  using  the  full 
width  at  half  maximum  (EWHM)  of  a  Gaussian  fit  to 
the  EDEs  in  Eigs.  8  and  9.  This  procedure  appears 
to  provide  a  reasonable  and  consistent  measure  of  the 
breadth  of  the  energy  distribution  and  avoids  the  un¬ 
certainties  that  would  otherwise  be  introduced  by  cal¬ 
culating  the  mean  and  standard  deviation  of  each  pro¬ 
file  in  Eigs.  8  and  9.  As  seen  in  Eig.  10,  the  widths  of 
the  ion  energy  distributions  are  measurably  larger  for 
the  low  magnetic  case  than  the  two  other  more  quies¬ 
cent  cases  with  higher  magnetic  fields.  Interestingly, 
the  highest  median  velocity  at  the  most  downstream 
location  was  for  the  must  oscillatory  case,  by  10  eV 
and  15  eV  relative  to  the  mid  and  high  magnetic  field 
cases.  Unfortunately,  it  is  not  precisely  known  which 
of  these  cases  produces  the  highest  efficiency.  Such 
measurements  would  determine  which  of  these  cases 
is  the  most  efficient  ion  accelerator.  At  present,  it  is 
surmised  that  the  more  oscillatory  cases  are  less  ef¬ 
ficient  since  more  high  energy  ions  are  created  deep 
within  the  acceleration  channel  and  are  likely  to  re¬ 
sult  in  more  ion  recombination  due  to  ion  wall  impact 
events. 

The  widths  of  the  profiles  are  better  compared  be¬ 
tween  the  three  cases  in  Eig.  11.  The  widths  of  the 
energy  distributions  for  the  nominal  and  mid  range 
radial  magnetic  field  strength  cases  are  nearly  the 
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Figure  10.  Median  ion  energies  bracketed  by  ion  distribu¬ 
tion  widths  for  the  three  magnetic  field  cases. 


same  through  out  the  acceleration  channel.  In  both 
cases,  the  widths  of  the  ion  energy  distribution  be¬ 
gin  near  zero  and  slowly  rise  reaching  a  maximum 
between  Z  =  -2  mm  or  Z  =  -1  mm.  Downstream, 
the  energy  distribution  widths  of  both  cases  appear 
to  remain  constant  at  approximately  40  eV  without 
appreciable  changes. 

Meanwhile,  the  behavior  of  the  low  (50%  of  nomi¬ 
nal)  radial  magnetic  field  strength  case  is  very  differ¬ 
ent.  The  energy  distribution  width  follows  the  other 
cases  until  the  Z  =  -4  mm  position.  Examining  Eig.  8, 
this  is  the  position  at  which  the  ions  first  experience 
any  significant  acceleration.  A  portion  of  the  ions 
have  already  been  accelerated  to  as  high  as  100  eV, 
and  an  even  larger  portion  have  just  been  ionized  and 
have  less  than  30  eV  of  acceleration.  In  contrast,  the 
less  oscillatory  modes  appear  to  have  already  expe¬ 
rienced  the  majority  of  their  ionization  and  gained 
approximately  50  eV  of  acceleration.  Downstream  of 
Z  =  -4  mm,  the  low  radial  magnetic  field  case  has 
a  rapidly  growing  energy  distribution  width  which 
grows  to  a  peak  of  175  eV  at  Z  =  -2  mm,  which  drops 
to  100  eV  at  the  exit  plane,  appears  to  drop  further 


8  of  10 


American  Institute  of  Aeronautics  and  Astronautics 
Approved  for  public  release;  distribution  is  unlimited. 


200eV  - 

+  Low  Magnetic  Field 

+  X  Mid  Magnetic  Field 

150  - 

"*"  +  O  Fligh  Magnetic  Field 

+  + 

100  - 

+  -1-  + 

+  +  + 

50  - 

8  is  8 

0  - 

*  *  ^  ^ 

T  '  I  '  I  '  I  '  I  f 


-8  -4  0  4  8mm 

Z  Postion 

Figure  11.  Measured  energy  distribution  widths  for  the 
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outside  the  thruster  reaching  a  width  of  70  eV  at 
Z  =  10  mm. 

One  possible  interpretation  for  the  evolution  of  the 
energy  distribution  of  the  nominal  and  80%  radial 
magnetic  field  strength  cases  is  that  an  axial  breath¬ 
ing  mode  oscillation  disperses  the  ion  creation  region 
by  a  potential  approaching  the  measured  width  of 
the  energy  distributions.  Alternatively,  it  is  interest¬ 
ing  to  note  that  the  width  of  the  two  more  quiescent 
cases  is  close  to  the  commonly  accepted  estimate  that 
the  peak  electron  temperature  within  the  acceleration 
channel  which  is  approximately  one  tenth  of  the  dis¬ 
charge  potential. 

The  coherent  oscillations  of  the  low  radial  magnetic 
field  case  may  manifest  as  a  similarly  coherent  axial 
breathing  mode  ionization  wave.  It  appears  to  be¬ 
gin  at  Z  =  -4  mm  rapidly  growing  and  peaking  at 
Z  =  -2  mm,  subsequently  decaying  in  amplitude  by 
Z  =  2  mm.  In  fact  at  Z  =  -2  mm,  the  ion  energy 
distribution  of  this  case  strongly  suggests  two  over¬ 
lapping  energy  distributions;  one  centered  at  190  eV 
and  another  at  70  eV. 

There  has  been  some  optical  evidence  that  high 
amplitude  oscillations  raise  the  electron  temperature 
significantly  in  a  Hall  thruster  anode  discharge. It 
appears  unlikely  that  the  electron  temperature  is  on 
the  order  of  the  the  170  eV  energy  distribution  width. 
An  increased  electron  temperature  likely  only  con¬ 
tributes  peripherally  to  the  increase  to  the  ion  energy 
distribution  width. 

Conclusions 

This  effort  has  non- intrusively  measured  xenon  ion 
velocity  profiles  in  the  center  of  the  acceleration  chan¬ 
nel  of  a  medium  power  (600  W)  Hall  thruster  from 
the  near  anode  region  to  10  mm  (0.15  exit  diam¬ 
eters)  downstream  of  the  geometrical  exit.  These 
measurements  were  taken  for  three  conditions  which 
varied  in  radial  magnetic  field  strength  in  a  ratio  of 


1.00  :  0.85  :  0.50.  During  these  tests  all  other  con¬ 
trols  (discharge  potential,  propellent  flow  rate,  back¬ 
ground  pressure,  cathode  heater,  and  cathode  keeper 
current)  were  unchanged.  In  the  three  cases,  the  an¬ 
ode  discharge  current  varied  by  only  2%.  The  pri¬ 
mary  observed  difference  between  the  three  cases  is 
the  substantial  increase  in  the  10  kHz  to  60  kHz  an¬ 
ode  current  breathing  mode  as  the  radial  magnetic 
field  strength  is  reduced. 

From  the  ion  velocity  profiles,  ion  energy  distribu¬ 
tions  are  computed.  The  case  with  highest  oscillation 
(and  the  lowest  radial  magnetic  field  strength)  has 
substantially  broader  ion  energy  distributions.  The 
LIF  diagnostic  allows  for  the  nonintrusive  character¬ 
ization  of  the  internal  plasma  oscillations  of  this  Hall 
effect  thruster. 

In  this  study  we  have  characterized  the  discharge 
current  oscillations  to  be  breathing  mode  oscillations 
generally  located  at  a  frequency  of  40  kHz.  The 
strength  and  spectral  width  of  the  breathing  mode 
oscillation  are  linked  in  the  conditions  studied  and 
also  appear  to  be  inversely  dependent  on  the  radial 
magnetic  field  strength. 

As  a  result  of  a  strong,  coherent  breathing  mode, 
the  ion  energy  distributions  are  up  to  4  times  broader 
than  when  there  is  a  weak,  frequency  distributed 
breathing  mode.  One  important  consequence  of  this 
is  that  the  highly  oscillatory  case  has  significant  pop¬ 
ulations  of  high  energy  ions  much  deeper  within  the 
Hall  thruster  than  the  quiescent  cases.  The  highly  os¬ 
cillatory  case  also  appears  to  have  an  extended  ion¬ 
ization  zone  within  the  thruster  acceleration  chan¬ 
nel.  Both  these  issues  are  most  likely  due  to  a  strong 
and  coherent  breathing  mode  axially  traveling  ion¬ 
ization  wave  which  oscillates  within  the  Hall  thruster 
both  distributing  the  ionization  region  and  dithering 
the  local  acceleration  potential.  It  is  believed  that 
this  two  effects  best  explain  the  broadened  ion  en¬ 
ergy  distributions  that  result.  It  is  also  possible  that 
the  relatively  high  frequency  of  the  breathing  mode 
may  result  in  an  elevation  of  the  electron  temperature 
within  the  Hall  thruster  acceleration  channel.  How¬ 
ever,  it  is  surmised  that  any  possible  electron  tem¬ 
perature  rise  only  has  a  minor  effect  on  the  breadth 
of  the  ion  energy  distribution. 

Since  the  strength  of  discharge  current  oscillations 
below  100  kHz  has  been  seen  to  increase  with  test 
chamber  background  pressure,  the  results  of  this 
work  imply  that  understanding  Hall  thruster  mea¬ 
surements  at  elevated  background  pressures  is  signif¬ 
icantly  more  complex  than  the  simple  ingestion  of 
background  neutral  atoms. It  is  clear  that  further 
investigation  of  these  phenomena  are  required. 

One  important  aspect  of  the  breathing  mode  os¬ 
cillations  is  their  possible  effect  on  Hall  thruster  life- 
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time.  The  lifetime  of  a  Hall  thruster  is  governed  by 
the  erosion  of  the  dielectric  walls  of  the  acceleration 
channel.  It  appears  from  this  work  that  long  term 
thruster  operation  in  a  highly  oscillatory  regime  may 
reduce  thruster  lifetime  due  to  the  introduction  of 
large  populations  of  high  energy  ions  deep  within  the 
acceleration  channel  combined  with  a  dithering  elec¬ 
tric  field.  We  are  not  aware  of  any  study  directly 
linking  oscillations  to  shortened  lifetime,  but  the  re¬ 
sults  of  this  work  support  this  supposition. 
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